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All the classes of static massless scalar field models available currently in the Einstein 
theory of gravity necessarily contain a strong curvature naked singularity. We obtain here a 
family of solutions for static massless scalar fields coupled to gravity, which does not have 
any strong curvature singularity. This class of models contain a thin shell of singular matter, 
which has a physical interpretation. The central curvature singularity is, however, avoided 
which is common to all static massless scalar field spacetimes models known so far. Our 
result thus points out that the full class of solutions in this case may contain non-singular 
models, which is an intriguing possibility. 
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Spherically symmetric solutions of Einstein equations for static massless scalar field configura- 
tions have been investigated in considerable detail in past. Bergmann and Leipnik l| were among 
the first to construct spherically symmetric static solutions for a massless scalar field. They had, 
however, only a limited success due to an inappropriate choice of coordinates. Around the same 
time Buchdal Q] developed techniques to generate solutions for this system, and also Yilmaz 
and Szekeres [J] found some classes of solutions for the static massless scalar field configurations 
in general relativity. 

Subsequently, Wyman |j| systematically discussed these solutions in comoving coordinates 
where the energy momentum tensor is in diagonal form, and showed a general method to ob- 
tain solutions in the case when the scalar field was allowed to have no time dependence. This gave 
a unified method to obtain most of the solutions obtained earlier. Also, Xanthopoulos and Zannias 
[fi] gave a class of solutions for time independent scalar fields in arbitrary dimensions where the 

ds conformally coupled to gravity have 
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spacetime metric was static. Further, the static scalar fie 
been a subject of immense interest to many researchers 7], 
fields have also been investigated in settings more general as compared to spherical symmetry (see 
e.g. H Q). 

The main interest in these models has been mainly due to several interesting properties that 



these solutions exhibit, as was pointed out, for example, by the JNW solution 



la ]. Mainly these 
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properties refer to the nature of the spacetime singularity and the event horizons in these space- 
times. There are no trapped surfaces in the model and the singularity which is visible at r = 2m 
has interesting properties [141 ] . 

These features were generalized in 15[, leading to a result that for static massless scalar fields, 
the event horizon is always singular in asymptotically flat spacetimes, and that these results are 
not necessarily restricted to spherically symmetric models only. In this sense, a study of static 
massless scalar fields coupled to gravity provides some important insights into the global structure 
of these spacetimes, and also it gives us useful information on the nature of singularities and 
trapped surfaces. 

The point here is that, the vacuum spherically symmetric model is the Schwarzschild solution, 
which is a black hole with an event horizon covering the singularity. However, an introduction of 
a smallest scalar field in the model radically changes the causal properties of the model, making 
the horizon and trapped surfaces to disappear and the spacetime singularity is visible. It is thus a 
matter of interest to examine if this class of models admit any singularity free solutions, in order 
to decide if the presence of a non-vanishing scalar field always causes a naked singularity. While 
this issue is examined here, in the process we also find a new class of models for static massless 
scalar field, indicating the possibility of non-singular massless static scalar field models. 

One of the main features of the Wyman class of solutions, which is the currently available most 
general class of models in the case under consideration, is that there exists in these spacetimes 
a central singularity which is naked. In the present note, we report a class of solutions where 
there is no such strong curvature singularity. However, a C° matching is necessary to achieve 
this, and as a result there is a shell which has singular matter. Many examples of this type of 
matching of spacetimes are available in the literature. Usually, this type of singularity can be given 
a physical interpretation, unlike the strong curvature singularity, and hence it is not considered to 
be pathological as the former. This type of matching conditions were first introduced by Lanczos 
and Isra el [lq| . Later on this has been used by many other authors (for a nice review on the topic, 
see e.g. [17(). The thin shell formalism has also been used for static spacetimes [181 ] . In our case 
also, a C° matching is performed in a static case. This class of solutions presented here is different 
from the Wyman class of solutions and does not have a naked curvature singularity. 

We consider here a four-dimensional spacetime manifold which has spherical symmetry. The 
massless scalar field <j){x a ) on such a spacetime manifold (M,g a b) is described by the Lagrangian, 
C = —^(j)- a (j)-bg ab ■ The corresponding Euler-Lagrange equation is then given by, 4>-ab9 ab = 0, and 
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the energy-momentum tensor for the scalar field, as calculated from this Lagrangian, is given as 



Let us consider the massless scalar field which is a Type I matter field [191 ]. i.e., the energy- 
momentum tensor admits one timelike and three spacelike eigen vectors. At each point q 6 M, 
we can express the tensor T ab in terms of an orthonormal basis (Eo, Ei, IE2, E3), where Eo is a 
timelike eigenvector with the eigenvalue p and E a (a = 1,2,3) are three spacelike eigenvectors 
with eigenvalues p a . The eigenvalue p represents the energy density of the scalar field as measured 
by an observer whose world line at q has an unit tangent vector Eg, and the eigenvalues p a represent 
the principal pressures in three spacelike directions E a . 

We choose the spherically symmetric coordinates (t, r, 0,(f>) along the eigenvectors (Eo, E a ), such 
that the reference frame is comoving, as was chosen by [5| and also [6j. The general spherically 
symmetric metric can now be written as, 

ds 2 = e 2 ^dt 2 - e 2 ^dr 2 - R 2 (t, r)dtf, (2) 

where dCl 2 is the metric on a unit 2-sphere and we have used the two gauge freedoms of two variables, 
namely, t' = f(t,r) and r' = g(t,r), to make the gt r term in the metric and the Ti r component of 
the energy-momentum tensor of the matter field to vanish. Thus the energy-momentum tensor has 
a diagonal form. We note that we still have two scaling freedoms of one variable available, namely 
t — > f(t) and r — > g(r). We note that the variable R represents the physical radius. 

We have for spherical symmetry <j) = <f>{t,r), and from equation (pQ) we see that T\q = <p'(j) = 0. 
So we have necessarily (j>(t,r) = 4>{t) or (f>(t,r) = (j)(r), with the energy-momentum tensor being 
diagonal. For the metric ([2]), and using the following definitions, 

G(t,r) = e- 2 ^(R') 2 , H(t, r) = e~ 2v (R) 2 , (3) 

F = R(1-G + H), (4) 
we can write the independent Einstein equations for the spherically symmetric massless scalar field 

n 



(in the units 8vrG = c = 1) as below (see [20(, [2l|]), 



p= i!ir (5) 
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u'(p + P r ) = 2(Pg-P r )^-P;, (7) 

-2R' + R'^+R^ = 0, (8) 

In the above, the function F(t,r), also called the Misner-Sharp mass, has the interpretation of 
the mass for the matter field, in that it represents the total mass contained within the sphere of 
coordinate radius r at any given time t. As noted above, in the static case, the metric components 
g^s are functions of r only necessarily, but the scalar field <f> itself can be in general either r or 
t dependent. In the case when (f> = (fi(r), which we consider here, the components of the energy- 
momentum tensor are given by, 

T\ = -T; = Tl = T* = ie~ 2 V 2 (9) 

It follows that the equation of state in this case, which relates the scalar field energy density and 
pressures is thus given by p = P r = —Pg. As noted by Wyman [5], there can be a class of static 
spacetimes where <f> = <j){t) also. But we would consider here only the class for which eft = 4>(r), 
which describes many earlier known interesting solutions for static scalar field spacetimes. 

We shall now consider the static spacetimes, when 4> = 4>{r) and g^ u = g^ u ( r )- The Einstein 
equations given above then reduce to the following set of equations, 

*" = ( *' _ it _ "''^ (12) 

e -2 tR ,2 = 1 _ F (13) 

R 

In the above, the equation (|12p can be integrated once with respect to r to give 



g-0— v+a 

where a = const. Eliminating now <j>' from these equations gives, 

I e -2u+2a 



(14) 



2 fl- H ' = F ' (15) 
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2 m R 2 R R 2 { j 

e-^R 12 = 1 - £ (17) 

We note that there is still a freedom left to transform the radial coordinate r, and thus the number 
of unknown variables is reduced to three in the three equations above. This freedom is just a 
coordinate transformation of the form 

r-> X (r), (18) 

which is allowed by the spherical symmetry of the spacetime. 

Apart from the above equations, we can also obtain a useful first integral, which is actually 
contained in the previous equations and it is advantageous to use it. The Einstein equations in 
this case can be written in the form R^ v = 0,„ <p, u . This implies i?oo = 0, from which we get 

. he*'" , s 

" - -sr (19) 

where the quantity h is a constant. From (fT9l) and (fill) , we get, 

= % + * (20) 

where a = he~ a and C\ are constants. 

To examine the Einstein equations above, we now define a function f(R) as below, 

f( R U=vp ( 21 ) 

The above is a general definition, and not any assumption, because the metric functions here 
depend on r only. Using this in (fTol) . we get, 

F = e 2a f(R) + C 2 (22) 

We can choose C 2 = 0, which gives F = e 2a f(R). Using this in (|17p . we get 



p -2i, d/2 _ i _ e a f( R ) 
R 



(23) 



Also, from ([2T]) we get, 



2u > = (Lm + 1)R> (24) 
J,R R 
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Taking e 2a = 1 and C2 = for simplicity and clarity of presentation, the last two equations, 
together with ([TBI) give 

R(f,R ? = (f - R) (f,R +Rf,RR ) " Rf,R (25) 

The equation above holds true in generality and we have not made here any assumption of a special 
coordinate condition or no specific radial coordinate choice has been made. The above represents 
clearly the main Einstein equation in the static case when cj> = (j)(r). The solutions to the same 
give the classes of allowed static massless scalar fields models in general relativity. The above is a 
non-linear ordinary differential equation of second order which is in general difficult to solve fully. 

We now consider a particular solution of (|25p . which is given by f(R) = — ^, which solves the 
above as is easy to check by inspection. This solution of f(R) gives a class of solutions to the 
Einstein equations as will be shown now. It should be noted here that for f(R) = — we have 
from d2IJ), e" 2ly = 2 and so v = const. Then from (I19p . we get h = 0, and therefore a = in (|20p . 
This shows that <j> can still be a non-constant function of r even when v is constant. 

From now on, we focus here only the class given by the condition f(R) = —1/R. To write down 
the solution explicitly and to specify the same in terms of the metric components, a choice of the 
radial coordinate r has to be made now. We note here that in the above consideration, there exists 
a freedom of choosing <f)(r), and therefore a choice for the same corresponds to actually making a 
choice of the radial coordinate r. Only after we have made a choice of (j>(r), would that scaling be 
fixed. 

In general, for the case / = — j^, from earlier equations we have, 

= ^V(r) (26) 



Then, the Einstein equation e 2 ^R' 2 = 1 — £ implies that we have, 

R 2 (R 2 + 1). 



R 

2R' 2 



(j)' 2 {r) 

By solving the above equation we get, 

R= ■ hf ,\. f ,s (27) 
smh(±-^0(r)) 

It thus follows that we can now write down explicitly all the metric components, in terms of the 
function 4>{r) and <f>'(r) above, thus giving a full solution. The metric coefficient v can be found 
out from (|2ip . The metric component R 2 can be found out from (I27j) . Finally the function ip is 
also known from (|26p . This completes the solution. 
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The function (j>(r) here can be viewed as a free function that we are actually allowed choose, and 
the energy conditions are always respected for all such choices. As an example, putting <j)(r) = -, 
we can recover a particular solution from the Wyman class Q|, which is given by, 

ds 2 = V - [ , J",' X . ]V - [ . - £ d Vdtf (28) 
2 smn (^i — i J sinn(-^i — i j 

We note that the range of (f> is from zero to infinity in this case. 

In general, we can in fact recover the entire Wyman class by this procedure. Each solution of 
(|25p with <f)(r) = - gives one solution from the Wyman class. The solution set of (|25p with this 
choice of 4>{r) gives the whole Wyman class of solutions. We note here that, in general, the choice 
of the function (j){r) is not just a gauge choice. For example, if we take 0(r) = — -, the range of <fi 
has an upper bound, but no lower bound and it goes from zero to negative infinity and the upper 
bound can be changed from zero to any other number also. In the choice that Wyman made, the 
range of 4> is different, and is necessarily restricted from below in that it goes from zero to positive 
infinity. Thus, in general the solution would be different from the Wyman class of models. 

For the metric given by (|28p . the Ricci scalar R c is given by 

Rc = - ( 29 ) 

Therefore there is a curvature singularity at the center R = 0. In what follows, we construct 
a class of solutions where the central singularity is absent, by choosing a specific form of (f>(r). 
Before proceeding further, however, we need to calculate the quantity Pi = f e^dr. This would be 
necessary to find out the proper length between two shells on any t = constant hypersurface. 

Pl(r) = [ e^dr = (30) 

Our purpose here is to construct and find a solution without any strong curvature singularity 
in the spacetime. Towards that purpose, we need to remove the central singularity R = 0, which 
for example exists in all other models, as discussed above. One way to achieve this is to construct 
a solution such that the physical radius R in fact does not vanish. To do this, we first notice, from 
(|27h . that R — > when <p(r) — > ±00. Also from (I30p . we see that Pi diverges when (j) —> 0. Since 
we are considering a static solution without any strong curvature singularity here, we can consider 
any t = constant hypersurface where the physical radius R must not vanish. In that case, there 
is no strong curvature central singularity in the spacetime. Further, since the spacetime is static 
here, it will be inextendible provided the proper length from any point on the spacelike hyperspace 
(any t = const surface) to the outer boundary of the hypersurface is infinite. This ensures the 
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regularity of the solution. If both of these requirements are to be satisfied, then <p cannot diverge 
and must go to zero twice in a range of the radial coordinate r. This implies that 4>{r) must have 
an extrema somewhere in that range of r where (j)'(r) = 0. 

There are in fact an infinite number of functions which satisfy these criteria, as required above, 
and so all of them are equivalent in this regard. Therefore, we make a simple choice here as given 
by, 

4>(r) = {r-a){b-r) (31) 

For this choice, (f> has a maxima at r = ( a ^ b ) where (j)'(r) = 0. From (I26p . it follows that = 
at this point, which means that the comoving coordinate system breaks down there. However, in 
the limit of r — > ^ & ^ , the curvature R c remains finite. Also, the proper length between the shells 
r = and r = b is infinite. So it is seen that this coordinate system covers the entire domain, 

Rmin <R< 00, (32) 

where R m i 



sinh 



l (b- a y 



l V2 4 

It is clear that the spacetime can be extended through the hypersurface R = R m in- We do this 
by joining two identical domains R m in < R < oo together at the hypersurface R = Rmin- In this 
case, there is no central curvature singularity in the spacetime and two such identical domains are 
glued together to give the full spacetime. To examine the matching at this joint, we need to find 
out the extrinsic curvature at the hypersurface r = . First, we rescale the time coordinate 
t' — V t, such that dr 2 = ^dt' 2 . We consider the orthonormal frame given by 

(33) 
(34) 
(35) 
(36) 

The extrinsic curvature of the hypersurface is then given by 

K (A){B) = -^ {A f V (B) V = n » e (A)^ W-B)' ( 37 ) 

with A,B = t, 9, 4>. 
In this case, 

| e 2^|coth(^=(r-a)(6-r)) 
K (6W) = -^p -R2 ' ( 38 ) 



?v = 


(0,|e 2 ^|,0,0) 




(1,0,0,0), 


e {B) ~ 


(0,0,l/i?,0), 


%) ~ 


(0,0,0, 1/R sin 
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and 



| e a0|coth(-^(r-a)(6-r)) 
K iv){v) = -#m> & ' ( 39 ) 



and the other components vanish. Now we note that, 



2-01 i 



lim 1 n , 1 coth(^=(r — a)(b — r)) = =F coth rJ- 



21 



(40) 



and thus the extrinsic curvature of the hypersurface r = r m i n is multi- valued. This means that 
there is a distributional source on this hypersurface, which we now calculate. The stress energy 
tensor of this source is given by Tu U = S uu S(rj), where rj is the Gaussian normal coordinate. 

S(A)(B) = ^ A f(B) = i( K (A)(B) ~ \ A)(B) trK + ) - {K- A){B) - h {A)(B) trK~) (41) 

where lim^ > ( a+b ) ±Q K(A)(B) = K(A)(B) an< ^ ^{A){B) is the induced metric on the hypersurface R = 

Rmin ■ 

From this it follows that, 

(b ct)^ 

S( T )(r) = -2 coth = -S {em = S m4>) (42) 

and the other components are zero. 

We note here that the surface energy density of this hypersurface is negative. We can, however, 
see now that the strong curvature singularity is removed in this case by introducing a C° matching 
at this hypersurface. The solution considered here is topologically different from the Wyman class 
of solutions, as there is no central shell in this case. There is a thin shell of singular matter in 
this case with a negative energy density. However, if we consider a sufficiently thick shell which 
includes the singular thin shell, then the mass inside that thick shell is positive. The negative 
energy thin shell then reduces the total mass content of the thick shell. In fact, this allows one to 
give a physical interpretation to the C° matching. While there is a mild singularity as discussed 
above at the joining surface, it can be given a sound physical interpretation (unlike the strong 
curvature singularity), as already discussed and used in the literature (see e.g. [22|). 

We also note here the fact that the spacetime we construct here is asymptotically flat. This can 
be seen in the following way. We consider only one part of the two identical domains, corresponding 
to Rmin < R < oo. Restricting ourselves only to this one part of the spacetime, we can consider the 
coordinate transformation so that the scalar field 4>(f) = i. This is the gauge that Wyman used to 
write down a class of solutions. However, our solution is different from that solution, because the 
coordinate system used by Wyman covers only a part of the spacetime manifold of the solution 
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given by us. Also, our solution does not have any strong curvature central singularity, unlike the 
Wyman model. In the Wyman coordinate system, the metric of our solution takes the form given in 
(|28p . From this expression, it is seen that in the limit f — >• oo, R — > oo, and then the metric becomes 
Minkowskian. So the part of the spacetime corresponding to R m in < R < oo, is asymptotically 
flat. Similarly, the other part of the spacetime can also be shown to be asymptotically flat. 
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